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ABSTRACT 

We analyzed near-infrared data of the nearby galaxy IC5063 taken with the Very Large Telescope 
SINFONI instrument. IC5063 is an elliptical galaxy that has a radio jet nearly aligned with the major 
axis of a gas disk in its center. The data reveal multiple signatures of molecular and atomic gas that 
has been kinematically distorted by the passage of the jet plasma or cocoon within an area of ~1 
kpc^. Concrete evidence that the interaction of the jet with the gas causes the gas to accelerate comes 
from the detection of outflows in four different regions along the jet trail: near the two radio lobes, 
between the radio emission tip and the optical narrow-line-region cone, and at a region with diffuse 
17.8 GHz emission midway between the nucleus and the north radio lobe. The outflow in the latter 
region is biconical, centered 240 pc away from the nucleus, and oriented perpendicularly to the jet 
trail. The diffuse emission that is observed as a result of the gas entrainment or scattering unfolds 
around the trail and away from the nucleus with increasing velocity. It overall extends for >700 pc 
parallel and perpendicular to the trail. Near the outflow starting points, the gas has a velocity excess 
of 600kms“^ to 1200kms“^ with respect to ordered motions, as seen in [Fell], Paa, or H2 lines. 

High H2 (1-0) S(3)/S(l) flux ratios indicate non-thermal excitation of gas in the diffuse outflow. 

Subject headings: Galaxies: active — Galaxies: nuclei — Galaxies: jets — Galaxies: ISM — Infrared: 

galaxies — Infrared:general — ISM: jets and outflows — ISM: kinematics and 
dynamics — ISM: lines and bands 


1. INTRODUCTION 

Gosmological simulations that required the implemen¬ 
tation oUWRnkdiolefe^bac^_to_j^gi^te_galaxy growth 
(e.g., iGroton et al.l 1200^ iBower et al.l l2006ll attributed 
a crucial role to radio jets. Jets can affect the inter¬ 
stellar medium (ISM) of galaxies by expelling, heating, 
and potentially compressing the gas. Heating was mainly 
invoked in the cosmological simulations because it pre¬ 
vents the accretion of new baryons from the cosmic web 
into galaxies, and because it disperses already accreted 
baryons into a hot, tenuous m edium, in which cooling 
is inefficient (|Maio et al.ll2007l ). Heating has been obser- 
vationally confirmed for the tenuous gas around Perseus 
A, from rip ples in X-ray-base d temperature and pres¬ 
sure maps ([Fabian et al.l l2006fl . It has also been pro¬ 
posed for the molecular gas in the outflow of 4C12.50, 
where the outflowing-to-ambie nt H2 mass is >30 times 
higher at 400 K th an at 25 K ([Dasvra fc Gornb^ 120121 : 
IDasvra et al.ll20i^ . Thermally driven bubbles drive gas 
expansion and, thus, produce the same observable signa¬ 
tures as shocks and ram pressure: outflows. 

Several outflow detection experiments have been suc¬ 
cessfully carried out for radio galaxies. Absorption from 
accelerate d H I gas at the t ip of radio hot spots was re¬ 
ported by iMorganti et al.l ([20051 1. Accelerated ionized 
gas close to the rad io jet lobes was detected in opti¬ 
cal w avelengths fe.g.. lRupke fc Veilleuxll2011t iShih et al.l 
I2OI3II . X-ray data revealed the presence of nuclear 
ultra-fast outflows, moving a t about 10% of the light- 
speed, in radio-loud galaxies ([Tombesi et al.ll2?)l^ . Re¬ 


cently, massive outflows of molecular gas were detected 
in galaxies of different morphological types showing 
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tections were considered of particular importance be¬ 
cause the molecular gas can probe star formation and 
because it can significantly contribute to the mass load 
of the outflow thanks to its high density. In some of these 
systems, the outflow entrained a few percent of the entire 
available reservoir. 

From a theoretical viewpoint, the detection of such 
massive outflows can be explained through the frag¬ 
mentation and reformation of dense molecular clouds 
because of the passage of tenuous ionized gas around 
them. Once the tenuous gas flow is initiated, thermo¬ 
dynamic (Rayleigh-Taylor and Kelvin-Hehnholtz) insta¬ 
bilities act upon the surface of the molecular clouds, 
dispersing them, an d facilitating the i r acceleration by 
ram pressure (e.g., iKlein et al.l 11994 iHookins fc ElvisI 
120101: l^ibovas fc Kindl2014l) . Molecules reform down the 
mass-loaded flow thanks to cooling, which happens on 
very short timesca les with respect to the flow travel time 
([Maio et al.ll2007ll . 

Following these results and their interpretation, mod¬ 
els treat ing radio jets as a heatin g source of the tenuous 
gas (e.g, iSiiacki fc SDringell2nn6[l were followed by mod¬ 
els treating radio jets as a source of s hock-induced energy 
dissipation upon the den se gas fe.g.. iWagner fc Bicknelll 
I2OIII : iGaibler et al.ll2012il . A major conclusion from the 
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work of I Wagner fc Bicknelll (l2011fl is that the jet cocoon 
influences the disk not only along the jet propagation 
axis, but also 500 pc (or more) in the transverse direction. 
This prediction needs to be tested in galaxies in which 
we can spatially resolve jets carving their way through 
disks. Such examples are rare according to the stan- 
dard unification schem e of active galactic nuclei (AGN; 
lUrrv fc Padovanil 1199511. One galaxy with this peculiar¬ 
ity is NGC 4258jCecil et al.ll200rf iWilson et ahll^OOll : 
ikrause et al.ll2007ti . IWilson et al.l (j2001tl associated dif¬ 

fuse X-ray emission 170 pc above the disk with hot gas 
in the jet cocoon. 

Another valuable example is that of 1C5063, a nearby 
elliptical galaxy with a central gas disk and a radio 
jet that propagates for severa l hundred parsec befo re it 
scrapes out of the disk plane (jMorganti et al.l[l998[ l. In 
radio data, bright synchrotron emission is seen at its two 
lobes, near the bases of the narrow- line-region (NLR) 
ionization cones (iMorganti et al.iri998[l . At the two radio 
lobes, as well as at a clump that is near the north-west 
edge of the NLR ionization cone, iKulkarni et al.l (Il998f l 
detected strong line emission in Hubble Space Telescope 
NICMOS narrow-band images that they attributed to 
gas excited by jet-related shocks. Strong polarization 
of the near-infrared (NI R) continuum light was seen a t 
the north radio lobe bv iLopez- Rodriguez et al. 1 (|2013ll . 
A fast neutral gas outflow was detected through H I 
absorption in front o f this lobe (IMorganti et al.l [19981 
lOosterloo et al.l 12000(1 . A molecular outflow, which is 
associated or thought to be ass ociated with th e same 
region, has been reported by IMorganti et al.l (l2013f l 
for the cold gas seen in millimeter CO data, and by 
iTadhunter et al.l (|2014fl for warm gas seen in NIR H 2 
data. These authors proposed that the outflow is driven 
by the impact of jet plasma upon dense ISM clouds. 

We analyzed SlNFONl data aiming to study the effects 
of the jet propagation on the atomic and, mainly, on the 
molecular ISM of 1C5063. We adopted Ho=70 kms“^ 
Mpc“^, 11 m= 0.3, and r2A=0.7 throughout our work. 

2. DATA REDUCTION 

The NIR data used in this work were taken from the 
archive of the European Southern Observatory (ESO). 
We reduced all available spectrophotometric data of 
1C5063 that were taken with the SlNFONl integral-field 
unit of the Very Large Telescope (VLT). On-source obser¬ 
vations of two hours, divided into 12 exposures of 600 s, 
were carried out in the H band on June 6 and June 7 2005 
(as part of the program 075.B-0348 A). Another two- 
hour-long observing run took place on May 22 and June 6 
2012, using the Hj-K instrument setup. This time, 24 on- 
source exposures of 300 s were taken (as part of the pro¬ 
gram 089.B-0971 A). All datasets were taken with a field 
of view (FOV) of ko"x8.0". The observations, which 
were carried out under seeing conditions of 1.0"±0.3", 
did not benefit from (laser-guide-star-assisted) adaptive- 
optic corrections. 

We started the data reduction by removing data rows 
affected by erroneous detector-hardcoded processing, us¬ 
ing the IDL routine provided by ESO for this purpose. 
We then used the ESO Recipe Execution (EsoRex) ver¬ 
sion 3.10 tool to generate a nonlinear-pixel map, a mas¬ 
ter dark frame, a master flatfield frame, and a bad- 
pixel map per night, using the routines sinfo_rec_detlin. 


sinfo_rec_mdark, and sinfo_rec_mflat, respectively. We 
then computed the optical distortion of the spectra of 
the multiple slitlets and the wavelength calibration solu¬ 
tion with the aid of lamp exposures (using the routines 
sinfo_rec_distortion and sinfo_rec_wavecal). We ran the 
last pipeline routine, sinfo_rec_jitter, to reduce individ¬ 
ual science frames of stars used as telluric/flux calibra¬ 
tors and of IC5063. This routine applied the previous 
calibration results to the science frames, and removed 
sky frames from the science frames. For the 2012 pro¬ 
gram data, one sky frame was assigned to two science 
frames because of the employed object-sky-object nod¬ 
ding strategy. 

Once the science cubes were reconstructed, an at¬ 
mospheric absorption correction was applied to them. 
To perform this correction, we used the spectra of 
HD 220787, HD 955, Hip063541, and Hip081166 for the 
nights of June 6 2005, June 7 2005, May 22 2012, and 
June 6 2012, respectively. Prior to scaling and remov¬ 
ing telluric spectra from science spectra, we identified 
and removed absorption lines that were intrinsic to the 
stars. Stellar cubes were then used for the flux calibra¬ 
tion. Namely, we opted to use HD 146332 and HD 220787 
for June 6 2005, HD 193933 for June 7 2005, Hip063541, 
Hip095404, and Hip098312 for May 22 2012, Hip081166 
and Hip096923 for June 6 2012. 

Our next task was to combine the flux-calibrated cubes 
of the different jitter positions and nights into a single 
mosaic. For this purpose, we aligned all images to the 
world coordinate system using the jitter offsets, occa¬ 
sionally correcting for small offsets pertinent to telescope 
pointing instabilities. We merged all available spectra, 
averaging flux values at any given wavelength and po¬ 
sition on the sky. Median instead of average flux val¬ 
ues were ascribed to pixels if the average value exceeded 
the corresponding median value by more than five times 
the standard deviation of the distribution. To obtain an 
ultra-deep (i.e., four-hour-long) exposure at the center 
of the mosaic, we combined all data after bringing the 
resolution of the H band cubes to that of the H-l-K band 
cubes (i?=1640). We trimmed outer parts of the mosaic 
with exposure times ;$1.5 hour, so that the signal-to- 
noise (S/N) ratios of the different regions are compara¬ 
ble. The final, fully reduced cube has a total FOV of 
8.3"x8.3" (or ^2kpcx2kpc). 

We identified spectral lines in the fully reduced cube, 
and removed the (stellar and AGN) continuum from each 
line within a range of 6000 km s“^. The continuum was 
calculated from neighboring spectral pixels along each 
spatial pixel. 

3. RESULTS 

3.1. Basic galaxy properties and regular gas motions 

The list of emission lines that we detected in the 
NIR spectrum of IG5063 comprises Paa, Bry, He I 
2.0587 /im, [Fell] 1.6440 /im, [Si VI] 1.9634 /im, [Si XI] 
1.9196 fim, [CaVIII] 2.321 /rm, H 2 (1-0) S(l), and H 2 
(1-0) S(3). A kinematic analysis was feasible for these 
lines, which were bright enough to be little affected by 
neighboring sky/telluric lines. The H 2 (1-0) S(3) line 
could be contaminated by He I and/or [Fe II] 1.954 ^m, 
—550kms“^ away (|Vanzi et al.ll200§l . Because we can 
neither test this assumption nor evaluate the strength of 
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these lines, we opt to show the H 2 (1-0) S(3) data, but 
to not use them in our kinematic analysis for velocities 
below —500kms“^. We detected numerous other emis¬ 
sion lines in the spectra, but we did not create kinematic 
maps for them, either because their flux was low or be¬ 
cause their profiles were substantially contaminated by 
sky/telluric lines. These lines comprise BrJ, BrlO, Bril, 
[Fell] 1.5330/rm, H 2 (1-0) S(0), S(2), S(5), S(7), and 
H 2 (2-1) S(l), S(3). The possibility that potential Brl2 
1.6412/rm emission significantly affects the [Fe II] 1.6440 
fim emission (even though the Brl3 and Brl4 lines are 
undetected) is evaluated and rejected in Section 15^ The 
H 2 ( 1 - 0 ) S(4) line was so contaminated by sky lines that 
its flux could not be recovered. Finally, we detected H 2 
(1-0) Q(l), Q(2), and Q(3) emission, but the correspond¬ 
ing lines were discarded from our kinematic analysis be¬ 
cause their restframe frequencies are only ~ 1000 kms“^ 
away and their fluxes are comparable: emission of am¬ 
bient gas in one line could be mistaken for emission of 
outflowing gas in another line. In the rest of this work, 
when referring to H 2 lines, we will refer only to lines 
probing transitions from the first to the fundamental vi¬ 
brational state unless otherwise noted. 

The emission line maxima at the location of the nucleus 
indicate that IC5063 is at a redshift of 0.0116±0.0006. 
Collapsed images of the H-band continuum and of 
lines within ±500 kms“^ from the systemic veloc¬ 
ity are shown in Fig [TJ The contours of these im- 
ages, as well as the contours of a 17.8 GHz image 
(|Morganti et al.l[2007[l . are plotted over an optical (547M 
filter) image take n with WFPC2 onb oard the Hubble 
Space Telescope (iSchmitt et al.l 120031) . The emission 
is spatially unresolved in [Si XI] and [CaVIII], which 
probe the galaxy nucleus. The stellar CO(3-l) ab¬ 
sorption seen near [Ca VIII] may be affecting the lat¬ 
ter’s distribution. In all other lines, the emission is re¬ 
solved and follows the radio jet trail, like the [O III] 
5007A emission in the WFPC2 image. [Fe II] shows 
the overall most diffuse emission. Four distinct regions 
can be identified in the NIR images: the nucleus (at 
20:52:02.40—57:04:07.68), the south radio lobe (1.8" or 
450 pc south east of nucleus, at 20:52:02.60—57:04:08.68), 
the north radio lobe (2.2" or 520 pc north west of the 
nucleus, at 20:52:02.17—57:04:06.55), and the extended 
clump near the north-west edge of the NLR ioniza¬ 
tion cone (3.8" or 900 pc away from the nucleus, at 
20:52:02.06—57 : 04:05 .05) that was previously seen by 
iKulkarni et al.l (|1998D . 

Spectra of these four regions are shown in Fig. |21 inte¬ 
grated over a radius of 0.5". The spectrum of the nucleus 
has a rising (from 2.4 mJy to 5.6 mJy) NIR continuum, 
characteristic of dust heated by the AGN. The highest 
number of detected lines is found in the nuclear region. 
Paa is the brightest of all, despite the strong atmospheric 
absorption window it is in. [Ca VIII], [Si VI], and [Fe II] 
1.6440 /im are among other bright lines. Contrarily, the 
region with the lowest number of detected lines is the 
NLR clump 3.8" from the nucleus. Its spectrum has a 
much flatter and fainter continuum, of ~0.4 mJy, as a 
result of the lower AGN and stellar emission. The slope 
and the intensity of the continuum at the north radio lobe 
and at the south radio lobe are intermediate of those at 
the nucleus and at the NLR clump. However, the H 2 


and [Fe II] emission is enhanced at the radio lobes - in 
particular at the north radio lobe, where these lines are 
brighter than at any other position in the central 2 x 2 
kpc^ of IC5063. 

The flux ratio of H 2 (1-0) S(l) over Paa is 2-8 times 
lower at the nucleus than within 100 pc from the north 
radio lobe, and <15 times lower further south of this re¬ 
gion (Fig. [3]). Likewise, this ratio is up to 2 times lower at 
the nucleus than within 100 pc from the south radio lobe, 
and <10 times lower further north of this region. The 
regions with emission excess from molecular gas, which 
could be star-forming or associated w ith shocked media 
([Ogle et al.ll2010l : iGuillard et al.l[2M^ . do not necessar¬ 
ily show emission excess from iron ions. The [Fell]/Paa 
flux ratio peaks mid-way between the radio lobes and the 
regions with H 2 emission excess, whereas the [Fe II]/ H 2 
(1-0) S(l) and [Fe II]/H 2 (1-0) S(3) flux ratios are max¬ 
imum in regions where the H 2 (1-0) S(l)/Paa flux ratio 
is minimum. High fractions of ionized to molecular gas 
emission are seen in regions elongated perpendicularly 
to the jet propagation axis. Such regions are likely to 
contain shocked gas, because Fe requires the passage of 
shock fronts to be liberated from dust grains and mas¬ 
sively deposited into a diffuse gas reservoir. The borders 
of the NLR are prominently outlined by ridges of high 
[Fe II]/Paa flux ratio, either due to shocks or AGN illu¬ 
mination. The clump near the NLR north-west edge is 
bright in [Si VI]. The logarithm of the [Si VI]/Br 7 flux 
ratio in that region is in the range 0 . 1 - 0 .3, which renders 
it compatible with shock and shock precurso r ionization 
models, but also with AGN photoionization (lAllen et al.l 
I 2 OO 8 I: [&oves &: Allen 11201^ iBicknell et al.ir2013l) . 

The effects of obscuration are evaluated by means of 
comparison of the Paa/ Bry ratio with the theoretically 
predicted value of 12 for case B recombination at lO^K. 
In most regions, the ratio takes values between 12 and 
9 (Fig. 131). These numbers respectively translate to op¬ 
tical depths of 0.0 to 2.2 in the H band, and 0.0 to 1.6 
in the K band, assuming that the absorption follows the 
Milky Way Rv=3.1 model of iDraini (I2003D . Data re¬ 
duction artefacts and differences in the gas excitation 
conditions can explain values outside this range. We 
find that there is little or no obscuration near the north 
radio lobe, near the nucleus, near the NLR most illumi¬ 
nated edge, and in front of the dust lane (north of the jet 
trail). Dust could be affecting the observed gas proper¬ 
ties ~300 pc west of the nucleus and near the south radio 
lobe, but any such effects are not quantifiable. Regions 
with high obscuration and H 2 / Paa excess coincide in 
proj ection with large-scal e structures, rich in molecular 
gas (IMorganti et al.ll2015D . 

To determine the kinematic properties of the gas in 
the center of IG5063, we fit a Gaussian function along 
the profile of each line in each spatial pixel. Only pix¬ 
els with S/N>10 were used. The results of this com¬ 
putation are presented in Fig. (H The disk previou sly 
known from HI observations (IMorganti et al.l I1998D is 
clearly seen in all of our tracers. It has a projected minor 
axis of 240(±30)pc, calculated from the He I map and 
from the hydrogen recombination line maps in regions 
displaying circular orbits. Its south-east side is moving 
away from us and its north-west side is coming toward 
us. On the south-east side, the average circular veloc- 
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ity Vc is 185±21kms“^ at 750 pc from the nucleus, and 
250±12 km s“^ at 900 pc from the nucleus. On the north¬ 
west side, the disk is rotating with —157±42kms“^, 
— 190±44kms“^, and —247±7kms“^ at 800pc, Ikpc, 
and 1.25 kpc from the nucleus, respectively. A confirma¬ 
tion of these numbers comes from the radio lobe spec¬ 
tra, where the jet-related shocks enhance the emission 
of the gas rotating with the disk. There, the shocked 
gas emits so prominently that disk-related spectral line 
displacements of order 200 kms“^ are visible in the scale 
of the entire NIR spectrum (Fig. [2]). The resolution- 
corrected velocity dispersion of the gas, as measured from 
the same gas tracers as those used for the rotational ve¬ 
locity measurement, is 115±20kms“^ at the nucleus and 
95±9kms“^ one kiloparsec north west of it. 

The examination of data cube slices around lines of in¬ 
terest reveals further details on the kinematics of the dif¬ 
ferent gas components. The disk rotation can be clearly 
seen in Fig. [SJ where each panel shows an image with 
a velocity bin equal to two spectral pixels of the initial 
cube: the bulk of the emission is shifted from the west 
of the disk to the east of the disk with increasing ve¬ 
locity. The emission from the north-west branch of the 
NLR peaks at about —320kms“^, and it can be seen for 
another ±300kms“^. Its counterpart on the other side 
of the galaxy, the south-east branch of the NLR, can be 
barely detected in the 157kms“^ and 315kms“^ Paa 
bins. Faint structures almost perpendicular to the disk 
major axis (at position angles between 25° and 50°) can 
also be seen at zero velocity in [Fe II] and Paa. One 
of them starts at the center of the galaxy and continues 
south west, whereas the other starts mid-way between 
the nucleus and the south radio lobe and continues north 
east. They are both linked to t he lar ge-scale spiral arms 
seen in CO bv iMorganti et all (l2015h . 

3.2. Gas near the radio lobes and along the jet trail: 
proof of its acceleration by the plasma 

A striking result in Fig. [5] is the detection of high- 
velocity gas near the projected positions of the radio 
lobes, where the jet is thought to interact with dense ISM 
clouds. Near the north radio lobe, the high-velocity flow 
is seen for multiple species (H, He, Fe, Si, H 2 ) and for all 
gas phases (i.e., atomic ionized/neutral and molecular). 
[Fe II], for example, is detected over a continuous veloc¬ 
ity range of 2000 km s“^ (Fig. [5]). The [Fe II]-emitting 
gas at ^1000kms“^ (Fig. [5]) is moving against the disk 
rotation by ~1200kms“^. In the direction of the disk 
rotation, about half of the [Fe II] emission comes from 
gas that is moving faster than the circular velocity along 
the line of sight (Fig. |2l bottom). Despite the potential 
existence of Brl2 1.6412/rm emission —510kms“^ away, 
we confirmed that the detection of the highly blueshifted 
[Fe II] wing is significant: under the assumption that the 
Brl2 flux can be as high as the Bril flux, we removed 
the stacked Bril image at —300±180kms“^ from the 
stacked [Fell] image at —810±180kms“^, and verified 
that the detection remains significant (Fig. [5]). Gas that 
is moving faster than 14 by more than 600kms“^ is seen 
in several other lines, including H 2 (1-0) S(l). 

Near the south radio lobe, the jet effects on the gas 
kinematics are only a little less pronounced than near 
the north radio lobe, even though the gas emission is 


lower by a factor of 2 (Fig. [51 lower panels). With 14 
rising quickly from ~100kms“^ to ~200kms“^ near the 
south lobe (Fig. |4|), we again detect gas moving with ve¬ 
locities 600kms“^ above the disk circular velocity and 
900kms“^ against it (Fig. [5]). This result, found for two 
sites, unambiguously shows that the irregular gas kine¬ 
matics are due to jet-cloud collisions - not due to some 
coinc idence (e.g., debris of a past merger; iColina et al.l 
B9911 in the line of sight behind one of the lobes). 

To study the outflow extent between the radio lobes, 
we created maps of the outflowing-to-ambient gas emis¬ 
sion (in [Fell] and H 2 ). We did so for two velocity 
ranges of the outflowing gas. With the fast outflow 
used in Fig. |6l we associated gas moving faster than 
600kms“^ (i.e., with a velocity exceeding I 4 by more 
than three times the disk velocity dispersion). With the 
intermediate-velocity outflow used in Fig. (3 we associ¬ 
ated gas moving slower than 600kms“^, but against the 
disk rotation. To be conservative, we used the maximum 
measured circular velocity, 250kms“^, and velocity dis¬ 
persion, 115kms“^ in our definitions. By ambient, we 
refer to gas moving slower than 14 (in either direction). 
Fig. ini shows that the fast outflow extends along the ra¬ 
dio jet trail from the vicinity of the nucleus to the north 
radio lobe. The intermediat e-velocity outflow ex tends 
further away (Fig. 0 see also IMorganti et ^120151 for a 
200kms“^ outflow of cold gas probed by CO). 

Fig. ini further reveals that the highly blueshifted 
and the highly redshifted emission near the north 
radio lobe are not cospatial. The emission at 
810(±)180kms“^ peaks 240pc south east of the emis¬ 
sion at—810(±)180kms“^ (Fig.Ej). 

To identify discrete outflow starting points, we used 
Hubble Space Telescope FOG optical images of high spa¬ 
tial resolution (~10pc). The FOG image shown in Fig. [5] 
is the average of three individual exposures probing po¬ 
larized light at 0 , 60, and 120 degrees, respectively. 
Taken with the F502M filter, these images comprise H/3 
and JG III] emission. Regions of interest are marked on 
Fig. [ 8 | With R1 and R2, we mark the beginning and 
the end of a 120 pc-long filamentary structure, which is 
perpendicular to the jet trail. The peak of the highly red- 
shifted [Fe II] emission is near R2. The dark region west 
of R2 is where the H-band optical depth is ^2, where the 
H 2 /Paa ratio is high (Fig. [3]), and whe re a large-scale 
CO s tructure reaches the inner gas disk (IMorganti et al.l 
HnH). With R3, we mark the NLR base, i.e., the region 
where the ionization cone becomes visible. With R4, we 
mark the peak of the highly blueshifted [Fe II] and H 2 
emission. There, the emitting region in the FOG data 
is oriented north-south (i.e., with an inclination with re¬ 
spect to the jet trail). Rl, R3, and R4, are respectively 
located 240, 320, and 400 pc away from the nucleus. 

Upstream the NLR cone and north west of R4, more 
outflowing gas is detected by its [Fe II] emission at 
315±90kms“^ (or ^550kms“^ against 14; Fig. [7|). A 
discrete starting point withi n this flow can be see n by 
the NLR clump identified bv iKulkarni et al.l (Il998f l. We 
mark it with R5 in Fig. | 8 j 

3.3. Diffuse gas out of dynamical equilibrium 

The flow of accelerated gas is not just parallel to 
the jet trail. Fig. [ 6 ] shows a diffuse outflow compo- 
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nent that is oriented perpendicularly to the jet trail, 
and that spreads south west of Rl. Its extent in ei¬ 
ther [Fe II] or H 2 exceeds 300 pc at high velocities and 
reaches ~500pc at intermediate velocities (Figs. [H [7]). 
The intermediate-velocity maps further reveal a similar 
diffuse outflow component spreading ~400 pc north east 
ofRl (Figs. [3®. The overall outflow has, thus, a bicon- 
ical shape. Both its cones extend beyond the projected 
disk thickness (as derived from regions with circular or¬ 
bits in the He I and hydrogen recombination line maps). 
Conical shapes of gas moving against the disk rotation 
can also be seen in the [Fe H 2 (1-0) S(l), S(3), and 
Bry kinematic maps of Fig. 0 

The center of this biconical outflow is located 240 pc 
north west of the galaxy nucleus. Its axis is perpen¬ 
dicular to the jet trail and parallel to the R1-R2 line 
in the FOC image, revealing an association between the 
outflow and the bright H/3 and [O III] emission. For 
[Fe II], the fraction of accelerated-to-ambient gas emis¬ 
sion exceeds one near this line (Figs. lU [Z])- It is worth 
noting that this is not necessarily the highest intrin¬ 
sic fraction, but the highest observed fraction due to 
ions scattered in our line of sight. The elongated re¬ 
gion with the highest observed accelerated-to-ambient 
gas emission continues to the north as a ridge of high 
[Fell] velocity dispersion (Fig. 0]). Along that ridge, 
the gas dispersion is boosted from 140(±10)kms“^ to 
290(±20)kms“^. On its west, the gas line-of-sight ve¬ 
locity goes down to —350(±20) kms“^ in [Fell], Paa, 
and [Si VI]. On its east, the iron ions are typically mov¬ 
ing with 180-300 km s“^. 

The data cubes provide enough information as to 
how the spatial distribution of the diffuse gas emission 
changes with velocity (Fig.jQj). In this Figure, the [Fe II], 
H 2 (1-0) S(l) and S(3) cubes are plotted at their finest 
velocity step, and the nuclear emission has been sub¬ 
tracted for visualization purposes. The diffuse emission 
from gas moving against the disk rotation is also visible 
prior to the nuclear flux removal (Fig. [5]). The inspec¬ 
tion of the individual cube channels reveals that the gas 
emission spreads from the nucleus to the north lobe with 
increasing velocity. The region that the diffuse gas occu¬ 
pies unfolds almost symmetrically with respect to the jet 
trail. At velocities >240kms“^, however, the symmetry 
breaks as the emission drops sharply north of the trail, 
near the ridge of high [Fe II] velocity dispersion. South 
of the trail, the emission continues, forming an arc that 
connects (in projection) the biconical outflow base with 
the north radio lobe. The arc is most pronounced for 
the molecular gas. An explanation for this progression 
in the outflowing gas spatial distribution as a function of 
velocity is given in the following Section. 

Position-velocity diagrams along position angles per¬ 
tinent to the biconical outflow are shown in Figs. ITOl 
and CII In Fig. Iini we examine a slice that is nearly 
perpendicular to the jet axis and that passes from the 
starting point of the biconical outflow. All along 800 pc, 
gas is continuously detected with a velocity peak of ^ 
200kms“^ for all of [Fe II], H 2 (1-0) S(l) and S(3). For 
[Fe II], we see that the emission spreads out to 1.2kpc, 
and that the gas attains a sharp terminal velocity of 
400kms“^. When reading this position-angle diagram 
from the jet propagation axis (i.e., from an offset of 3") 


to outer regions (i.e., to offsets of 0" or 6"), we find that 
the diffuse gas velocity increases with increasing distance. 
The velocity gradient is smooth, and close to symmetric 
in the two sides of the jet. 

In Fig. [11] we examine slices nearly parallel to the jet 
axis. Along the northern border of the biconical out¬ 
flow, which runs parallel to the jet propagation axis at 
a distance of ^350 pc, there is intermittent [Fe II] emis¬ 
sion at 900kms“^ along both sides of the disk iFisf. [TT|) . 
This indicates that the ionized gas wind generated north 
west of the nucleus is strong enough to dominate the 
[Fe III emission even north east of the nucleus (see also 
Fig-Sj). This is not true for the warm molecular gas, even 
though H 2 (1-0) S(3) emission is also seen along 1.3 kpc 
in both sides of the disk (at about — 100kms“^). West 
of the nucleus, part of this emission is coming from the 
disk. East of the nucleus, the emission is blueshifted by 
260kms“^ with respect to the disk rotational velocity. It 
is, thus, likely due to shocks generated east of the nucleus 
(Figs, m ITT]) . The region with negative velocities north 
east of the nucleus can also be seen in the H 2 (1-0) S(l) 
and S(3) velocity fields (Fig. |4|). This is another region 
where gas flows in both line of sight directions, as [Fe II] 
emission is also detected there at ^900kms“^. Diffuse 
outflowing gas is, thus, also present on the east side of 
the disk. Its emission is, however, too weak to reveal 
changes in the gas spatial distribution as a function of 
velocity. 

Notably, the (more tenuous) ionized gas attains 
higher velocities than the (denser) molecular gas at 
fixed positions (Figs. E] [9] (tO] [TT]) . This result is in 


agre ement with multi - phase outflow acceleration mod¬ 
els (iHookins fc Elvis 120101 iWagner fc Bicknelll 120111 : 
iZubovas fc Kindl2014h . 


4. DISCUSSION 


4.1. On outflow properties pertinent to galaxy evolution 

A question pertinent to galaxy evolution studies that 
NIR data help us to address is whether the jet heats up 
the molecular gas it accelerates. To evaluate this possi¬ 
bility, we compared the gas excitation temperature Tex 
in the ambient medium and in the intermediate-velocity 
diffuse outflovQ. For this purpose, we created temper¬ 
ature maps of the warm H 2 using the (1-0) S(l) and 
S(3) lines, assuming that the gas can be described by 
a single temperature component in local thermodynamic 
equilibrium (LTE). The results are shown in Fig.fT^ The 
ambient gas excitation temp erature at the north r adio 
lobe is 1500±200K (see also iTadhunter et al.l 120141) . In 
the rest of the disk, it is 2000±200K. At the edges of 
the disk, near the regions with excess of molecular-to- 
atomic-gas emission (Fig. [3]), it goes up to 2700±500K. 
The low temperature of the ambient gas at the north 
radio lobe could be indicative of a region that is denser 
and more self-shielded than its surroundings, which can 
be described by several components of subthermally ex¬ 
cited gas. 

For the diffuse accelerated gas at 310±70kms“^, Tex 
varies significantly from one region of the map to an¬ 
other. In the south of the jet trail, Tex is within the 


^ This question cannot be evaluated for the fast outflow due to 
the potential blending of the H 2 (1-0) S(3) line with [Fe II] or He I 
lines at —550kms“^. 
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range of values found for the ambient medium. How¬ 
ever, no Tex solution can be found for many pixels north 
of the jet trail: north of the regions R3, R4, and the 
north radio lobe (where diffuse 17.8 GHz emission ex¬ 
ists; Fig. [5]), the (1-0) S(3)/S(l) flux ratio is above 
the saturation limit (2.04) for a population that is de¬ 
scribed by a Boltzmann distribution and that decays 
spontaneously (Fig. [12]). This result implies that the 
molecular gas is in part non-thermally excited and that 
radiative decay in the form of fluorescence is plausi¬ 
ble. The non-thermal source of excitation could be X- 
rays or UV light coming, for example, from the jet- 
related high-velocity shock wave or from the AGN ra¬ 
diation (iLepp &: McCravlll983t [B l ack fc van Dischoeck I 
Il987t iDraine fc Woods! 119901 : iMouril 11993b This result 
does not rule out a temporary temperature increase of 
the gas that got exposed to the non-thermal source of 
excitation: radiative cascade and heating are often cou¬ 
pled processes (e.g., lMourilll994[). Evidenc e of heating 
comes from the work of lYoung etld] (|2007D . who found 
that the dust is warmer at the north radio lobe than at 
south radio lobe. 

The H2 mass distribution that we can create when tak¬ 
ing the Tex map at face value and assuming an ortho- 
to-para ratio of 3 is shown in Fig. jT^j The total am¬ 
bient gas mass is 3500 Mq (of which 990 Mq, 460 Mq, 
and 150 Mq are found within a radius of 0.5" from the 
north lobe, the nucleus, and the south lobe, respectively). 
For the gas at 310±70kms“^, the total mass is 170 Mq 
( of which 30 Mq and 80 Mq are within 0.5" from the 
north lobe and the nucleus, respectively). The remaining 
60 Mq come from the diffuse medium. Fig. |T3| demon¬ 
strates that the maximum mass of the diffuse molecular 
gas (that can be detected from our line of sight) fol¬ 
lows a filamentary distribution near the biconical out¬ 
flow axis. The filamentary structure carries 35 Mq. For 
comparison, the mass of the molecular gas in the fast 
outflow at —810(±)210kms“^ would be 25 Mq, if all of 
the H2 (1-0) S(l) emission in it was due to thermal- 
ized gas at Tex=2000K. The total mass in the outflow 
at all velocities and regions is, thus, u nlikely to exceed 
some hundreds of solar masses (see also iTadhunter et al.l 
1201411 . We, thus, conclude that the contrast between the 
mass carried in the warm phase and t he mass carried 
in the cold phase (2.3x 10^-1.3x 10® Mq; iMorsfanti et al.l 
I2OI3II remains stark. If the acceleration of the molecular 
gas happens through its partial dissipa t ion an d reforma¬ 
tion as proposed by iHookins fc Elvis! (1201011. then this 
resul t is a proof that microturbulence (iNesvadba et al.l 
l2011tl ceases rapidly enough for cooling and potentially 
for star formation to occur during the outflow propaga¬ 
tion timescale. 

4.2. The jet and its cocoon as drivers of the gas 
kinematics 

The bottom line of the results presented in Section O is 
that numerous gas components with irregular kinematics 
are detected in the central square kiloparsec of IC5063. 
Which mechanism drives the observed gas motions? In 
this Section, we argue that the propagation and the de¬ 
flection of the radio jet plasma does so. 

The diffuse gas geometric characteristics are not in line 
with an AGN-related or a starburst-related mechanism 
of acceleration. The biconical outflow in IG5063 cannot 


be associated with a bubble driven by the AGN radiation 
(as in NGC4258: iJimenez- Vicente et al.ll2010[l . because 
it is not centered at the nucleus. The gas kinematics are 
also incompatible with those of an expanding bubble: 
we observe no spherical symmetry in the velocity gradi- 
en t of the accelerated Fe ions. No ultra-fast outflow (as 
in iWagner et a.l.ll2013ll has been de tected in X-ray spec¬ 
tra of IG5063 ( Tombesi et al.ll20T^ . and no [Ga VIII] or 
[Si XI] outflow is seen in the SINFONI data. A nuclear 
starburst is ruled out for the same symmetry reasons 
as those applicable to the AGN. An obscured starburst 
located near region R1 of Fig.|8]is also unli kely. A bicon¬ 
ical outflow similar to th at in, e.g., M82 (iWalter et al.l 
I2OO2I: iVeilleux et al.ll200^ would mainly show redshifted 
emission in the northern hemisphere of the galaxy and 
blueshifted emission in the southern hemisphere of the 
galaxy. The SINFONI data indicate a symmetry of the 
redshifted emission in the north and south hemispheres, 
and an asymmetry of the redshifted vs. blueshifted emis¬ 
sion to the west and east of region R3. 

Contrarily to an AGN or starburst driven wind, nu¬ 
merous findings support the jet-ISM interaction: the gas 
emission is brighter north west than south east of the nu¬ 
cleus (like the plasma emission itself), the non-thermally 
excited gas lies parallel to the jet trail, outflows are ini¬ 
tiated in at least four discrete regions along the jet trail, 
the highest gas line-of-sight velocities (<—600kms“^ 
or >600kms“^) are attained near the proposed jet- 
ISM interaction points, and the gas emission can be 
blueshifted/redshifted in the same hemisphere. 

As the plasma propagates, bow shocks push the sur¬ 
rounding tenuous gas and form an expanding cocoon. On 
its turn, the cocoon applies ram pressure to clouds that 
it encounters as it expands. The emission from the lat¬ 
ter clouds can be either blueshifted or redshifted on the 
same side/hemisphere of the galaxy, depending on the 
positioning of the jet and the gas with respect to us. If 
the jet travels between us and background clouds, then 
we see redshifted emission from the background clouds 
that are being pushed further away from us by the co¬ 
coon. If foreground clouds exist between the jet and 
us, then we see blueshifted emission from the foreground 
clouds t hat are being pushed close r to us by the cocoon 
(see also lGarcia-Burillo et al.f2014fl . In regions where the 
jet plasma directly impacts clouds, then both blueshifted 
and redshifted emission from the scattered material can 
be detected. If the jet propagates through a clumpy disk, 
the plasma can get deflected on several dense regions and 
the cocoon can push several clouds. Then, the cloud dis¬ 
tribution is most important for determining how many 
outflows will b e initiated and in which directions they 
will be headed (IWagner fc Bickn^l201lD . 

One schematic representation of the jet-ISM geometry 
explaining all observed outflows in IC5063 is presented 
in Fig. [131 It follows the cloud naming convention of 
Fig.lU under the assumption that the large-scale H/3 and 
[O III] emission in the F OG image is associated with the 
jet trail, as advocated bv iMorganti et all (1200711 . Fig. [TT] 
shows a jet propagating through a clumpy medium and 
its cocoon on the disk plane. North west of the nucleus, 
most clouds are located behind the jet, whereas south 
east of the nucleus, most clouds are located in front of 
the jet (with respect to us). This is in line with several 
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results: with the small inc lination of the jet with respect 
to the disk plane (Fig. [8] iMorganti et al.l [200711 making 
the north west radio lobe visible below the dust lane in 
the southern galaxy hemisphere, with the NIR obscura¬ 
tion being higher at the south east radio lobe than at the 
north west radio lobe (Fig.|3]), and with the gas emission 
in the FOC image being diffuse south east of the nucleus 
- possibly due to high ram pressure (Fig. [8]). 

In this picture, the jet and its cocoon drive a mostly 
redshifted outflow north west of the nucleus and a mostly 
blueshifted outflow south east of the nucleus, clearing 
out some of the gas in the disk. North west of the nu¬ 
cleus, consecutive major jet-ISM interactions are occur¬ 
ring, starting at the base of the biconical outflow (near 
Rl; Fig. [T4l) . There, the passage of the plasma by the 
near side of a dense I SM structure, as that seen in the 
ALMA CO(2-l) data (jMorganti et al.lF2ni5h . causes the 
redshifted biconical emission. Part of the plasma is likely 
to be scattered to the south of the disk plane, as indicated 
by the weak 17.8 GHz emission, which extends along the 
southern biconical outflow axis and whic h ends where the 
fast outflow ends (jMorganti et al.l[2007[ l. Then, near the 
north radio lobe, in region R4, an interaction between 
the jet and the far side of another gas clump generates 
the outflow coming toward us (Fig. [TTl) . This clump, 
oriented north south, could have been shaped by plasma 
that got deflected at the north radio lobe and that drilled 
a channel on its wa y out of the disk (Fig. [8]) . Again, the 
17.8 GHz image of IMorganti et al.l (120071 1 shows diffuse 
radio emission north of region R4 (starting north of R3 
and ending near R5). Between regions Rl and R4, at 
the NLR base (R3), we see either the spatial overlap of 
the outflows at Rl and R4, or another jet-ISM interac¬ 
tion that scatters gas in many directions. Further up 
the NLR cone, in region R5, the jet again pushes away 
from us the near side of a cloud. On the other side of 
the nucleus, near the south radio lobe, scattered gas is 
accelerated in both lines of sight near the tip of the radio 
jet (Figs. IHl lMj) . 

The progression in the location of the diffuse gas emis¬ 
sion with velocity can be explained by the jet propagation 
in more than one way. We could be observing tenuous 
ISM material that got accumulated and accelerated be¬ 
fore impacting the dense structure at Rl, material that 
got scattered from that dense structure, or a combina¬ 
tion of the two. In the first case, we observe the motions 
of a cocoon-driven outflow that gets mass loaded as it 
moves from the nucleus to the outer parts of the galaxy. 
The velocity of this outflow will be higher near the jet 
cocoon (i.e., closer to the dense structure at the biconi¬ 
cal outflow starting point). Once at the dense structure, 
the mass-loaded outflow will appear deflected in two di¬ 
rections perpendicular to the initial trajectory because 
of line of sight effects. In the second case, we observe 
an outflow that got mainly mass loaded after the im¬ 
pact of the jet upon the dense structure, when gas got 
bounced to different directions. The gas scattered close 
to our line of sight appears to be moving faster and to 
be located closer to the scattering surface, whereas the 
gas scattered nearly perpendicular to our line of sight 
appears to be moving slower and to be closer to the nu¬ 
cleus. In either way, scattering of the plasma links the 
gas entrainment along the jet trail and perpendicular to 
the jet trail, as well as in small and large scales. The 


high-velocity-dispersion regions near Rl to R4 (e.g., the 
ridge north east of Rl in the [Fe II] kinematic map) are 
indicative of scattering. A spatial overlap of two accel¬ 
erated/scattered gas components in projection or in vol¬ 
ume (i.e., at the interface between two flows) can add to 
the observed dispersion. 

The diffuse outflowing gas that we detected within the 
central square kiloparsec of IG5063 reaches hundreds of 
parsec away from visible jet breaking points. Its total 
north-east to south-west extent exceeds 700 pc. In the 
south of the jet trail, the outflowing gas emission fills 
the entire volume between the galaxy minor axis and the 
NLR south-west branch (Fig. [5]). Likewise, in the north 
of the jet trail, the regions with an intermediate-velocity 
outflow (Figs, m |8]) together with the region of oversat¬ 
urated H 2 (1-0) S(3)/S(l) emission (Fig. [12]), again fill 
the volume between the galaxy minor axis and the NLR 
north-west branch. This emission, thus, spreads for 700- 
900 pc along its two south-east to north-west borders. 
The outflow detection in an area of order 1 kpc^ empha¬ 
sises the role of jets in affecting the r eservoir for star for¬ 
matio n. Notably, the simulations of I Wagner fc Bicknelll 
(|2011f l mimicking the propagation of a jet in a clumpy 
disk, predicted that the turbulence would be so wide¬ 
spread. 

The distance out to which we detect the diffuse gas 
depends on the flow travel time or on the presence of a 
physical obstacle such as star-forming regions (with ex¬ 
cess H 2 /Paa emission; Fig.|3|), in addition to the depth 
of the observations. The gas travel time gives us a limit 
on the timescales of the (intermittent or not) relativistic 
plasma ejection. If the Fe ions were moving at an aver¬ 
age projected velocity of 900kms“^, then 4x10® years 
would be required for them to travel from the disk to 
350 pc north of it. For the molecular gas moving at 
— 100kms“®, the timescale would be an order of mag¬ 
nitude higher, but it would be more uncertain: the frac¬ 
tion of molecular gas that forms in the flow from dense 
atomic gas and the pertinent timescales for its forma¬ 
tion are unknown, rendering the use of a constant, low 
velocity unclear. 

5. SUMMARY 

We used SINFONI data of IG5063, a nearby elliptical 
galaxy in which a radio jet propagates through a disk, to 
look for unambiguous observational evidence of the jet- 
ISM interaction and to study its imprints on the ISM. 
Multiple pieces of evidence were found, which can be 
summarized by the detection of high-velocity gas in (at 
least) four regions along/near the jet trail, and of diffuse 
accelerated gas that is associated with the jet cocoon. 
Main results are as follows. 

• Emission is seen all along the spatially resolved 
jet trail in numerous NIR lines, including molec¬ 
ular and atomic (neutral/ionized) gas tracers. The 
brightest H 2 and [Fe II] emission is located at the 
projected position of the north radio lobe. The 
south radio lobe and a clump in the NLR north¬ 
west branch are clearly identified. Diffuse atomic 
and molecular gas emission is also seen between the 
nucleus and the north radio lobe. 

• At the north radio lobe, we detect gas moving with 
>700kms“^ above the disk circular velocity and 
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1200 kmagainst it. Likewise, at the south radio 
lobe, we detect gas moving with 600kms“^ above 
the disk circular velocity and 900kms“^ against it. 
The derivation of this result for both sites unam¬ 
biguously proves that the gas has been accelerated 
by a collision of the jet plasma with ISM clouds. 

• We also detect an outflow bright in [Fe II] emis¬ 
sion at 315kms“^ that starts near the north west 
branch of the NLR, and another outflow bright in 
[Fe II] and H 2 emission that starts 240 pc north 
west of the nucleus and that is associated with a 
filamentary structure perpendicular to the large- 
scale jet trail. This structure, which is character¬ 
ized by strong H/3 and [O III] emission and sur¬ 
rounded by extended 17.8 GHz emission, could be 
due to deflected plasma. The newly detected out¬ 
flow near this structure has a biconical shape. It 
entrains iron ions and H 2 molecules moving as fast 
as —8I0±180kms“^ and 310±70kms“^, respec¬ 
tively. Along its southern axis, it shows the highest 
mass load in the line of sight, and higher [Fe II] 
emission than in the ambient medium. 

• The spatial distribution of the diffuse outflowing 
gas changes as a function of velocity from 0kms“^ 
to 480 km s“^. More specifically, the [Fell], H 2 (1- 
0) S(l) and S(3) emission unfolds parallel to the jet 
propagation axis and extends further away from the 
nucleus with increasing velocity. Likewise, the ve¬ 
locity of the diffuse gas smoothly increases with in¬ 
creasing distance from the biconical outflow start¬ 
ing point both to the north east and to the south 
west of the jet trail. Again, plasma propagation 
and scattering explain the observed gas kinematic 
behaviour. 

• Gas is driven out of dynamical equilibrium and 
dispersed over an area of order Ikpc^ by the 
jet and its associated cocoon. In the west side 
of the disk, outflowing gas is detected in either 
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phase (atomic/molecular) or direction for more 
than 700 pc along the jet propagation axis. Parallel 
to this axis, but north and south of it, the emission 
from gas that is kinematically or excitationally af¬ 
fected by the jet fills the entire volume between 
the galaxy minor axis and the NLR. This emission, 
thus, spreads for ^Ikpc. The total north-east to 
south-west extent of the diffuse outflowing gas also 
exceeds 700 pc. Diffuse gas exists in the east side 
of the disk too. Its extent cannot be determined, 
however, because the diffuse emission from acceler¬ 
ated gas is weak there, like the radio emission. 

• The highest line-of-sight velocities are attained 
near the areas where the plasma gets scattered. In 
the flow, the atomic gas typically reaches higher ve¬ 
locities than the molecular gas at fixed positions. 

• The computation of an excitation temperature 
from the H 2 (1-0) S(l) and S(3) line fluxes, under 
the assumption that the gas is thermalized, fails 
for gas in the diffuse outflow near the north radio 
lobe. This result implies that the gas is in part 
non-thermally excited and that radiative decay in 
the form of fluorescence is plausible. It doesn’t, 
nonetheless, rule out gas heating. 
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Fig. 1.— Left: images of the H-band continuum and of bright lines detected in the SINFONI data. The line images were created 
by collapsing continuum-free channels over the velocity range -500kms“^<V<500kms~^. Right: contours of these images (at 10a level 
steps going up to a m aximum of lOOcr) are plotted over a Hubble Space Telescope WFPC2 547M image comprising [O III] 5007A emission 
JSchmitt et alj[2003l) . To indicate the radio jet projected position, contours of a 17.8 GHz image IIMorganti et al.ll2007l) are also overplotted 
in the top right panel. 
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Fig. 2.— Top: rest-frame NIR spectra, extracted within a radius of 0.5^^ from the nucleus, the north radio lobe, the south radio lobe, 
and the bright NLR clump of IC5063. A wavelength-independent constant of 2 mjy, 1.5 mjy, and 1.5 mjy has been added to the spectra 
of the north lobe, of the south lobe and of the NLR clump, respectively, to facilitate their visual comparison. Bottom: continuum-free 
spectra of the same regions, zoomed in at spectral lines of interest. All panels are plotted using the same color scheme, i.e., black for the 
nucleus, blue for the north radio lobe, cyan for the NLR clump, and red for the south radio lobe. 
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Fig. 3. — Line ratio maps, created from the images presented in Fig. ^ (and comprising the emission from gas with velocities 
—500kms“^<V<500kms~^). Crosses mark the gas emission maxima at zero velocity at the nucleus and near the radio lobes. 
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Fig. 4.— Flux, velocity, and instrumental-broadening-corrected velocity dispersion of various gas tracers, as indicated by the best-fit 
solutions to the continuum-subtracted line spectra in each spatial pixel. Crosses mark the gas emission maxima at zero velocity at the 
nucleus and near the radio lobes. The dashed gray lines in the velocity field of [Fe II] outline the biconical outflow perpendicular to the jet 
propagation axis. 
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FIG. [4]— continued. 
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Fig. 5. — Continuum-free slices of the SINFONI data, binned by a factor of two. Empty frames correspond to channels that have been 
masked because of the presence of sky/telluric residuals or neighboring ISM lines that could be confused with high-velocity wings. Crosses 
mark the gas emission maxima at zero velocity at the nucleus and near the radio lobes. 
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FIG. [5]— continued. 
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Fig. 6.— Ratio of the outflowing gas emission over the ambient gas emission. The fraction is computed for velocity bins of lkms~^, 
because of the different velocity ranges spanned by the two components. Velocities in the range —250<V<250 kms“^ were associated with 
the ambient medium. For the gas in the outflow, we used velocities <—600kms“^or >600kms“^ (exceeding the disk circular velocity by 
more than three times the disk velocity dispersion). The outflow velocity range is noted in the upper left corner of each panel. Signal below 
three times the noise of the collapsed image has been masked. At negative velocities, the maximum predicted Brl2 contribution has been 
removed from the [Fe II] image. 
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Fig. 7.— As in Fig. 0 but now using the [Fell] emission at 315(di90)kms ^ (left panel) and the H 2 (1-0) S(l) emission at 
310(di70) kms~^ (right panel) instead of in the fast outflow. Regions east of the nucleus have been masked, as they can be contami¬ 
nated by the regular disk emission at these velocities. 



Fig. 8 . — Regions of interest with respect to the radio jet trail (left panel) and the outflow (middle and right panels) are plotted over a 
high angular resolution Hubble Sp ace Telescope FOC im age of IC5063. The FOC F502M filter was used, which comprises and [O III] 
emission. 17.8 GHz contours from [Morganti et al.l II2007I) are plotted in blue in the left panel. The north lobe and south lobe are marked 
with NL and SL. Contours of the outflowing-to-ambient [Fe II] emission (from Figs. (S) 0 are plotted in the middle panel, revealing 
regions with a high fraction of outflowing Fe ions. Cyan, red, and magenta curves correspond to highly blueshifted, highly redshifted, and 
intermediately redshifted emission, respectively. The contours are at steps of 0.1 starting from 0.1 for the fast outflow, and at steps of 
0.25 starting from 0.25 for the intermediate-velocity outflow. In the right panel, we plot a map of the outflow in either gas phase (atomic 
or molecular), using pixels where >10% of the emission comes from outflowing gas. The flow extent at intermediate velocities is shown 
in magenta. The fast outflow extent is shown in red for the region with gas flowing away from us, in blue for the region with gas flowing 
toward us, and in purple for the region with gas flowing in both directions. 
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Fig. 9. — Data frames showing the spatial progression of a faint structure from the nucleus toward the north radio lobe with increasing 
velocity. The disk emission, approximated by a Pao image at 0ili75kms“^ that was normalized to the peak nuclear emission, has been 
removed from each frame to assist the visualization of the diffuse gas. 
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Fig. 10.— Position-velocity diagram along a slice almost perpendicular to the jet propagation axis that passes from the [Fe II] biconical 
outflow starting point. The offset along the slice increases from north east to south west and the flux is normalized to unity. The absorption 
at negative velocities in the [Fe II] panel is real, due to the stellar CO(7-4) bandhead. The absorption at —400kms“^ in the H 2 (1-0) S(3) 
panel is due to a sky line. 



- 57 ® 04 ' 05 " 


- 57 '* 04 ' 10 " 


[Fell] (0+135) km/s 



20 ^ 52 '" 02?5 02?0 


1 


0 



[Fell] 



Fig. 11.— Position-velocity diagram along slices almost parallel to the jet propagation axis. The offset along the slices increases from 
south east to north west and the flux in each panel is normalized to unity. Letters are used to indicate which position-velocity diagrams 
correspond to which slices. 
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Fig. 12.— Maps of the warm molecular gas excitation temperature in the ambient medium (left panel), and in the range 240kms ^ 
to 380kms“^ comprising diffuse gas that moves against the disk (right panel). The maps were calculated under the assumption of gas in 
LTE. In the right panel, the excitation temperature of the gas in the dark red pixels either exceeded 5000K, or it could not be computed 
under the assumption of thermalized gas. Emission from the disk could be contributing to the flux of both lines. For this reason, the side 
of the disk moving away from us has been masked. 
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Fig. 13.— Maps of the warm molecular gas mass in the ambient medium (left panel), and in the range 240kms“^ to 380km s~^ 
comprising diffuse gas that moves against the disk (right panel). The maps are calculated for the excitation temperature maps of Fig. 1121 
In the right panel, emission from the disk could be contributing to the measured H 2 mass. 




observer 


Fig. 14.— Schematic representation of the jet-ISM interaction. Solid black lines represent the radio jet. The gray area depicts the cocoon 
in the plane of the disk. Clouds located behind the radio jet with respect to the observer are pushed further back by the jet plasma/cocoon. 
Their outflow is characterized by redshifted emission (red arrows) in either galactic hemisphere. Inversely, clouds located in front of the 
radio jet are pushed nearer to the observer (blue arrows). The region nomenclature follows that of Fig. [8] 


























